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has been carried out. The mean and the variance of the total kinetic energies of the fission fragments 
have been calculated and compared with the experimental values. 
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I. INTRODUCTION 

At present it is commonly agreed upon that the fission process is a dissipative phenomena, where initial energy of 
the collective variables get dissipated into the internal degrees of freedom of nuclear fluid giving rise to the increase in 
internal excitation energy causing evaporation of a large number of precision neutrons, which can not be explained by 
standard statistical model calculations. As dissipation is referred to the interaction of the system coordinate with the 
large number of degrees of freedom of the surrounding reservoir, this process is always associated with the fluctuation 
of relevant physical observable. This mesoscopic description is inevitable once the fluctuations of the observables are 
amenable to experimental observation. 

Keeping this fact in mind, the fission process is picturised as follows; the collective variable such as elongation axis 
acts as a ’Brownian particle’ interacting stochastically with large number of internal nucleonic degrees of freedom 
constituting the surrounding ’bath’. Several workers solve either the Langevin equation 01 , or multidimensional 
Fokker-Planck equation 00], which is a differential version of Langevin equation, in order to study the time evolution 
of probability distribution function and have calculated the mean as well as the variance of the total kinetic energies 
of the fission fragments. 

The experimental data show that the root mean square deviation of kinetic energy is always very small compared 
to the mean kinetic energy (~ 0.1). The question naturally arises whether one can utilise this simple fact in the 
theoretical scheme instead of solving the Langevin equation or corresponding Fokker-Planck equation in detail. In 
this spirit, we assume that the full solution of the Fokker-Planck equation admits an asymptotic expansion in terms 
of strength of the fluctuation. Admittedly, in the zeroth order of fluctuation one should obtain macroscopic behavior 
of the dynamical evolution. Obviously, when the fluctuation is ignored one obtains the deterministic picture and the 
machinery to handle such situation is to solve the Euler-Lagrange equation of motion. This was done in our previous 
works H 0], where dissipation was generated through nonconservative Rayleigh function and the total kinetic energies, 
the fission yields and neutron multiplicities were calculated. However, in the above scheme, there had been no scope 
to study the effect of fluctuations originating from the stochastic dynamics of fission on fission observables. Here, 
we report a modified scheme where fluctuations have been included in the fission dynamics by making asymptotic 
expansion of the probability distribution function in terms of intensity or strength of the fluctuations as argued 
before. Thus, the picture we adopt here, is as follows: Due to smallness of the relative fluctuation, the process is 
grossly described in terms of macroscopic equation. The stochastic description is introduced by studying the evolution 
of the probability distribution of narrow width over its deterministic values. 

In what follows, in Sec.II we describe briefly the procedure and derive the equations of the moments. The results 
of the calculations are discussed in Sec.III. Finally, concluding remarks are given in Sec.IV. 


II. ASYMPTOTIC EXPANSION OF THE FOKKER-PLANCK EQUATION 

A. The method 

As argued before, the mesoscopic description of the fission process begins with a set of Langevin equations: 

X i = h i ({X},{Y}) + r ]i (t) (2.1a) 

Yi = Hi({X}, {V}) ;*= 1,2,3, ...,N (2.1b) 

where hi and Hi are given functions of the stochastic collective variables X\, X 2 , ..., Xn and Yi, Y 2 ,..., Y)v in the 
fission process and refers to the driving noise term associated with the interaction of the «th collective variable 
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with the reservoir constituting nucleonic degrees of freedom. For simplicity, we assume the noise to be a gaussian 
white with zero mean and decoupled for different degrees of freedom with auto-correlation functions given by 


< r]i(t) >= 0, < r]i(t)rij(t l ) >= Di(yi)5(t - t')Sij, 


( 2 . 2 ) 


where £h(yj) is the diffusion coefficient associated with ith variable, depending only on the sample space yi for the 
stochastic variable Y, i: . 

The Fokker-Planck equation corresponding to the Langevin equation (2.1) is 


df{{x},{y},t) 

dt 


d{hjf) 
^ 1 dxi 


d(HJ) 

dyi 


(1/2)A(2/O0] 


(2.3) 


The quantity f({x},{y},t) is the probability density function depending on the variables aq, X 2 , ..., xjv, yi, j/ 2 > •••, Vn 
and time t explicitly. If we are interested in finding the time evolution of the conditional probability distribution 
function then we have to solve Eq.(2.3) with initial values x,(0) = x®, yi{ 0) = y?, Vi, at t = 0. That is, we have to 
solve Eq.(2.3) for those realisations which are known to start from these specific points in the whole sample space. 

The asymptotic expansion was developed by van Kampen || for constant diffusion constant. The method consists 
of writing the stochastic variables as the sum of deterministic value and a fluctuating part at each time t with root 
of the diffusion constant as a strength of the fluctuating part. In the present paper we generalise this method for the 
situations where the diffusion coefficients depend on the stochastic variables explicitly. Such a situation is encountered 
in the case of fission process, where the friction coefficient depends explicitly on the collective variable or shape of the 
nucleus at each instant of time. In this case, we further assume that, in the asymptotic expansion, the strengths of 
the fluctuating parts of the stochastic variables depend only on the deterministic values of the respective y variables: 


Xi = Xi + V D (yi)(i (2.4a) 

yi = Vi + V D (.yi)€i (2.4b) 

The quantities {Ci}, {6i} refer to the fluctuations of the stochastic variables {xi} and {y.;} around their deterministic 
values {xi},{yi}. Next, we introduce the new distribution function Q depending only on the variables {C*}>{£i} and 
t. The normalisation condition suggests that the / and Q will be related by 

N 

/(M, m, t) = n [A^r^aa m, t) (2.5) 

2 = 1 


Substituting Eq.(2.4) in the Fokker-Planck equation(2.3),making Taylor expansion of h({x}, {y}), H({x}, {y}) 
around {x},{y} and collecting coefficients of various order of D(yi ),we could generate a hierarchy of equations.As 
expected,the first set would give rise to the equation of motion for {x} and {y}. 

Xi = hi({x},{y}) (2.6a) 

yi = Hi({x},{y}) ; V* (2.6b) 


Eqs.(2.6) are the Euler-Lagrange equation for deterministic motion.These equations are to be solved with initial 
conditions {x(0)} = (x°},{y(0)} = {y 0 }. Next, we are going to calculate the conditional probability distribution 
f({x}, {y}, t | {a: 0 }, {y 0 }, 0) or Q({C}, {£}, t \ 0, 0,0). 

Assuming the variation of diffusion coefficient over the narrow width of the distribution function at any instant of 
time to be Q(D), we could replace the second Fokker-Planck coefficient D(y) by D(y) at each instant of time.This 
assumption makes the calculation extremely simple.Collecting coefficients O{D 0 ), we get back quasilinear Fokker- 
Planck equation for Q: 


dQ 

dt 


+£( 

i 


D{Vi) 

D{Vi) 


)Q 


d(c* Q) , , dfoQ) , d&Q) , J d(0Q) 

V I1, — +b, — +c ‘^r +di ~ar 


<1/2) f] 


where a,;, bi,Ci,di are given by 


, 9h _ D(yi) 
^xj [ 2 D(yt) 
dh 

w 


(2.7) 


(2.8a) 

(2.8b) 
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( 2 . 8 c) 


dH D( Vi ) 
« = (te-) - (7 


v 2£>(yi)' 


, f dH 

* = ( »> 


Eq.(2.7) suggests that 




( 2 . 8 d) 


(2.9) 


where the distribution function Qj for each j satisfies the similar equation written below without the subscript: 


3Q , ,*>(»),„ , a«0) , k 9({Q) , 8(S) , /((«) „„,3 2 0, 

(—)« = + t— + c— + d— - W2) w \ 


dt K D(y) 
subject to the initial condition 

The solution of Eq.(2.10) is given by 

Q((,Z,t) = r 1 


Q(c,z,t = o) = 6(cm) 


exp - {i*< + Ut + m# + GW* + 2c(t)kl] }dkdl 


l (2t r) 2J 7 J ^ ' " s ' 2 D(t) 

where g(t),G(t),c(t) satisfy the set of coupled first order differential equations : 

g dh dh D 

2 = <s> 9+( 5; )<:+ 2 

G <J!f ,, ,3H 

¥ = ( ^ )G + < & )c 

dh ,dh ,dH ctff 

c=( s )c+( ^ )G+( ay )<:+( & )9 

with the initial conditions 

5 (0) = G(0) = c(0) = 0 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13a) 

(2.13b) 

(2.13c) 

(2.14) 


Once Q(C)C, t) is known,from Eq.(2.9) and Eq.(2.5) the full conditional probability distribution function 
f{{x},{y},t | {x 0 }, { 5 0 }, 0) is known. Integrating this function over all variables except one,say Xi, one identifies 
5 i(f) as the variance of the stochastic variable X.;. 


< (X,- < X, >) 2 >= 9 l (t) 


(2.15) 


We note that the homogeniety of Eq.(2.10) suggests that < C,(t) >=< £(t) >= 0, or the average of the variables 
X and Y at any time will be determined by their deterministic values obtained by solving Euler-Lagrange equation 
(2.6). Similarly one observes from Eq.(2.12), 


< (X- < Xi >) >< (Yi- <Yi>)> = a(t ) 

< (Yi- < Yi >) 2 > = Gi(t) 


(2.16a) 

(2.16b) 


B. Application to the fission process 

In the fission process,we choose the elongation axis and velocity associated with it as the stochastic variables 
interacting with a large number of internal nucleonic degrees of freedom constituting a heat bath at temperature T 
determined by the excitation energy available to it. We further assume that the ’collisional’ time scale of the nucleonic 
degrees of freedom is much shorter than the time scale of the macroscopic evolution of the collective variable so that 
at each instant of time the heat bath is assumed to be in quasi-stationary equilibrium. 
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As stated before the Euler-Lagrange equation(2.6) was solved in our earlier works [||. To avoid repetition we 
deliberately omit the procedure and scheme to solve that equation. For the sake of completeness we merely write that 
equation and refer to our previous paper to clarify the details. 

Giving correspondence to the terminology used in this paper, we associate 


Y = r,X = 


(2.17) 


Thus we have 


H(x,y) = x = r 

u , , r^ 2 • d(Vc + V N ), f 

K r ’ *0 = “ T - - Q- r -1/m 




(2.18a) 

(2.18b) 


The quantities Vc,Vn represent the Coulomb and nuclear interaction potentials, respectively. The quantity 7 ,g 
and L refer to friction coefficient, reduced mass associated with the fissioning liquid drop and the relative angular 
momentum, respectively 

While solving for variances we appeal to Eq.(2.13). In that set the diffusion coefficient D is evaluated employing 
Einstein’s fluctuation dissipation theorem. The set thus reads 


r)h f)h 

m = 2 (-)g(t) + 2 (-)c(f) + 2 7 (r)T(r)/^ 2 
G(t) = 2 c(t) 

... ,dh. . . „,dh. . 

c (t) = (gv)c(i) + 2(—)G(t) + g(t ) 

with the initial conditions(2.14). The initial conditions of r and r for solving Eq.(2.18) are 

( E*R n , 


r(t = 0 ) = r min , f(t = 0 ) = (- 


2 fi 


1)1/2 


(2.19a) 

(2.19b) 

(2.19c) 

( 2 . 20 ) 


where the fissioning nucleus starts from the minimum of the potential energy surface, the quantity Rn is a random 
number between 0 and 1 from uniform probability distribution and E* is the available energy. 

Solving Eq.(2.18) and Eqs.(2.19) simultaneously with initial conditions (2.14) and (2.20) we generate the conditional 
probability distribution function f(r,r,t \ r(t = 0 ),r(t = 0),0). The probability distribution function f(r,r,t) could 
be obtained as 


f(r, r,t) = j f(r , r, t\r(t = 0 ), r(t = 0 ), 0 )f(r(t = 0 ), r(t = 0 ), 0 )dr(t = 0 )dr(t = 0 ) 


( 2 . 21 ) 


where f(r(t = 0 ),r(t = 0 ), 0 ) is the probability distribution of position and velocity of the stochastic variables at 
the initial time. As described by the initial condition(2.20), this can be represented as 


f(r(t = 0 ), r(t = 0 ), 0 ) = 8(r(t = 0 ) - r min ) x f(r(t = 0 )) 


( 2 . 22 ) 


Here, we assumed that each fissioning nucleus in the ensemble starts from a fixed initial position but with different 
partioning of initial excitation energy 0. Finally, substitution of Eq.(2.22) in Eq.(2.21) would give 


f{r, r,t) = ^2 f(r, r, t \ r min , ( E ) 1 / 2 , 0 ) 


Rn 


2 g 


(2.23) 


III. RESULTS AND DISCUSSIONS 


The temporal evolutions of the variables g(t), c(t), G(t) along the fission trajectory have been computed by solving 
numerically the set of Eqs. 2.19c. The results are plotted in Fig. [j] for a representative system 16 0 + 124 Sn. It is 
seen from the figure that, initially, all of them increase steeply and then their magnitudes become nearly constant 
throughout the rest of the trajectory. Furthermore, the calculation shows that 
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c 2 (t) 

g{t)G{t) 


< 1 


(3.1) 


This implies that the correlation of position and velocity of the elongation variable (r) is much smaller compared 
to their respective variances. This fact simplifies our calculation of variance of energy of the fission fragments. The 
variance of energy and average of kinetic energy at scission point are given by 

= (/*■ r) 2 g(t ) + 2 G(t) (3.2a) 

< E (t) > = g-gitsc)/ 2 + Edet (3.2b) 


where t sc is the time at scission point and Edet is the deterministic value of total fragment kinetic energy (TKE) 
after scission and ~ 100 - 200 MeV. In deriving the above approximate results Eqs. (3.2), we utilise our observation 
(3.1). It is further assumed that the variation of the potential over the narrow width of the probability distribution 
is small so that the average of the potential is approximated as the value of the potential at the mean position. The 
variation of the kinetic energy variance a E as a function of time has also been displayed in Fig. [I], The value of cj 2 e 
is also seen to increase steeply at the beginning and then it becomes nearly constant throughout the rest of the time. 
As envisaged earlier, the result clearly shows that ^ E> <C 1, which demonstrates the validity of asymptotic expansion 
in deriving the result instead of solving the Fokker-Planck equation in detail. 

The theoretical predictions of <jE{th) for the fission of several compound systems produced in the 158.8 MeV 18 0 
and 288 MeV 16 0 reactions on various targets have been displayed in Fig. ||alongwith the respective experimental data 
p| for comparison. The experimental data aEifixp) is represented by the filled circles and the solid curves are the 
results of the present calculati ons. It is seen from Fig. |] that when the projectile energy (and vis-a-vis the excitation 
energy of the fused composite) is relatively lower (lower half ), the calculated values are in fair agreement with the 
data. However, the calculation underpredicts the experimental value of <je for the heaviest target considered ( 238 U 
in the present case). With the increase in the projectile energy (and the excitation energy of the composite) (upper 
half of Fig. |2j), the theoretical predictions are found to underestimate the corresponding experimental values and the 
discrepancy between the two increases with the increase in mass number. We have also studied the fragment mass 
asymmetry dependence of energy variance, a E (A\, A 2 ) for some representative systems and the results are displayed 
in Fig. It is seen from Fig. || that the theoretical values of variances have only a weak dependence on the fragment 
mass asymmetry. 

In order to investgate into the discrepancy between the predicted and the experimental values of the TKE variance, 
it is observed that the experimental values of <jE(fixp) are usually obtained by averaging over the full mass yield 
spectrum. Therefore, aE(exp) consists of two terms, viz.. ( i ) contributions arising due to stochastic fluctuations in 
the dynamics of fission process and ( ii ) contributions from the variation of the mean kinetic energy with the fragment 
mass asymmetry. So, aE{exp) may be written as [ fTTI , 

erf; (exp) = cr| + a%(kin) (3.3a) 

4 = 5>|(Ai,A 2 ).y(A 1 ) (3.3b) 

Ai 

cr E (kin) = Y^IE- < E(Ai,A 2 ) >] 2 • Y(Ai) (3.3c) 

Ai 


where <j 2 e (Ai : A 2 ) and < E{A\, A 2 ) > are the variances and mean values of the total kinetic energy of two fission 
fragments with mass numbers A\ and A 2 (compound nucleus mass Acn = A\ + A 2 ), E being the average of < 
E(Ai,A 2 ) > over the normalised fragment mass yield, Y{A\) with 4 Y(A\) = 1 . It is, therefore, clear that for 
a proper comparison of the theoretical predictions of kinetic energy variance with the relevent experimental data, 
the theoretical numbers should be averaged over the respective fragment mass distribution. Moreover, the calculated 
value of cr E (th) should be compared with the extracted experimental value of a E obtained by substracting a E (kin) 
from a E (exp). 

We have calculated a E (kin) for a few systems for which the experimental fragment mass yield data are available 
@> taking < E(Ai, A 2 ) > from Viola systematics ]l^|. The results are given in Table | alongwith the values of a E 
extracted from the experimental data. It is evident from the Table | that the contribution to the variance from the 
variation of mean TKE with fragment mass asymmetry, a E (kin ), increases with the increase in excitatation energy 
and mass of the composite. For the lighter system e.g. O + Ag : it is observed that this contribution is quite large 
due to the asymmetric nature of fission. Further, it is seen from the Table Q that there is a very good agreement 
between the UE{th) and the <je extracted from experiment for the systems for which the mass fragment yield data 
are available. Moreover, the values <je, i.e., ( a E (exp) — cr E (kin ) ), are also shown in Fig.|| as open triangles and they 
agree very well with the predicted values of TKE variance. 
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IV. SUMMARY AND CONCLUSIONS 


We have developed asymptotic expansion of the Fokker-Planck equation for the systems in which the relative 
fluctuation of collective variable is small and have derived the equation of variances. The formalism is applied to the 
case of fission where the fluctuation in total kinetic energy is small as compared to its mean value. The calculated 
variances, aE(th ), underpredict the experimental data, aE(e.xp). The discrepancy between the two increases with the 
increase in the mass of the composite and the excitation energy of the composite. However, we emphasize that this 
discrepancy is not very surprising because the experimental values contain additional contributions from variation of 
the mean kinetic energy over full fragment mass yield distributions. Once this contribution is properly taken care of, 
the predicted TKE variances are found to agree quite well with the <je extracted from the experimental data for the 
systems where the fragment mass yield data are available. Therefore, for a more direct test of theoretical models it 
is necessary that experimental estimation of variances should not have admixture of other contributions arising due 
to the variation of mean kinetic energy over differnt mass yields. This may be achieved if measurements are done in 
smaller mass bins. 

Thus, it may be concluded that the present approach is quite successful in reproducing the extracted TKE variances 
from the experimental data without going in to solving the Fokker-Planck equation in detail. In the present studies, 
the correlation of the position and velocity of the elongation axis has been found to be small. However, in the cases 
where such condition is not valid the energy variance still can be calculated by adding a term 2pr( df ' Vc ^ VN ^ )c(t). 

The procedure developed here could systematically generate higher order hierarchies for relatively larger fluctuations 
than the ones encountered in the present studies. In those cases one may have to solve the higher order equation 
which would involve higher order derivatives of the functions h(x,y ) and H(x 7 y), in general. 
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FIG. 1. Variation of g(t), c(t), G(t) and a% as a function of time t for the system 16 0+ 124 Sn. 


FIG. 2. Variation of cte as a function of target mass number, for 288 MeV 16 0 (upper half), and 158.8 MeV ls O (lower 
half) induced fission reactions. Filled circles correspond to the experimental data (Ic| and solid lines correspond to the present 
theoretical results. Open triangles are the modified results using Eq. (3.3) (see text). 


FIG. 3. Variation of predicted values of <je as a function of fragment mass asymmetry, a syrn = |Ai — A 2 \/Acn- 


TABLE I. Calculation of TKE variance 


System 

E lab 

(MeV) 

a%(exp) 

(MeV 2 ) 

S Al (< E(A u A 2 ) > —E) 2 Y(Ai) a 
(MeV 2 ) 

(MeV 2 ) 

cT e 

(MeV) 

CJE(th) 

(MeV) 

18 O + iM Sm 

158.8 

112.3 

10.1 

102.3 

10.1 

9.7 

ls O + 197 An 

158.8 

190.4 

15.2 

175.2 

13.2 

11.9 

le O + 197 An 

288.0 

331.2 

49.7 

281.5 

16.7 

13.3 
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16 0 + 109 Ag 288.0 225.0 142.8 82.2 9.0 

a See Eq.(3.3) 
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